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ABSTRACT 

Context. The heavy mass loss observed in evolved AGB stars is usually attributed to a two-stage process: atmospheric levitation by 
pulsation-induced shock waves followed by radiative acceleration of dust grains, which transfer momentum to the surrounding gas 
through collisions. In order for an outflow to occur the two stages of the mass-loss scheme have to connect, i.e. the radiative accel- 
eration can only be initiated if the levitated gas reaches a distance from the stellar photosphere where dust particles can condense. 
This levitation distance is limited by the kinetic energy transferred to the gas by the shock waves, which imposes strict constraints on 
potential wind-driving dust species. 

Aims. This work is part of an ongoing effort aiming at identifying the actual wind-drivers among the dust species observed in circum- 
stellar envelopes. In particular, we focus on the interplay between a strong stellar radiation field and the dust formation process. 
Methods. To identify critical properties of potential wind-driving dust species we use detailed radiation-hydrodynamical models 
which include a parameterized dust description, complemented by simple analytical estimates to help with the physical interpretation 
of the numerical results. The adopted dust description is constructed to mimic different chemical and optical dust properties in order 
to systematically study the effects of a realistic radiation field on the second stage of the mass loss mechanism. 
Results. We see distinct trends in which combinations of optical and chemical dust properties are needed to trigger an outflow. Dust 
species with a low condensation temperature and a NIR absorption coefficient that decreases strongly with wavelength will not con- 
dense close enough to the stellar surface to be considered as potential wind-drivers. 

Conclusions. Our models confirm that metallic iron and Fe-bearing silicates are not viable as wind-drivers due to their near-infrared 
optical properties and resulting large condensation distances. TiC>2 is also excluded as a wind-driver due to the low abundance of Ti. 
Other species, such a Si02 and AI2O3, are less clear-cut cases due to uncertainties in the optical and chemical data and further work 
is needed. A strong candidate is Mg 2 Si04 with grain sizes of 0.1-1 yum, where scattering contributes significantly to the radiative 
acceleration, as suggested by earlier theoretical work and supported by recent observations. 
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1. Introduction Observations tell us that there exists a number of different 

dust species in the circumstellar environment (see, e.g.,[Molster 
|et al.| |2010| for an overview), however, most of these are not 
abundant enough or have insufficient cross-sections to generate 
the radiative acceleration needed to drive a wind. Furthermore, 
in order for a dust species to trigger a wind it has to form close 
to the stellar surface, within reach of the shock waves. Due to 
an environment dominated by a strong radiation field the closest 
distance from the star at which a given dust species can exist 
will depend on both the optical and chemical properties of the 
material. 

The chemical composition of dust grains in the circum- 
stellar environments of AGB stars is first and foremost de- 
termined by the elemental abundances in the atmosphere. As 
stars evolve during the AGB phase their atmospheric chemistry 
may transform from oxygen-dominated (M-type, C/O < 1) to 
carbon-rich (C-type, C/O > 1) due to ongoing nuclear burn- 
ing and dredging-up of processed material (see, e.g., the review 
by [Herwig 2005 ). Since most of the carbon and oxygen in the 
extended atmosphere will be tied up in the tightly bound CO- 
molecule, whatever element is more abundant will dominate the 
dust chemistry. As a consequence, carbon-based condensates are 



The prevailing scenario for stellar winds observed in AGB stars 
is a two-stage process, involving atmospheric levitation by shock 
waves and radiative acceleration of dust grains. In the first step 
pulsation-induced shock waves catapult material outwards in al- 
most ballistic trajectories. The pulsations themselves are not suf- 
ficient to cause an outflow that is comparable with observed 
mass loss rates or velocities (e.g. [Wood! |1979| |Bowen| [1988 ) 
but the shocks typically manage to lift the gas to a few stellar 
radii. As a consequence of the propagating shock waves in the at- 
mosphere, density-enhanced layers of gas form around the star; 
layers where the temperature is low enough for dust formation to 
occur (e.g Fleischer et al. 1992 Jeong et al. 2003 1. The dust is 
of crucial importance in the second step of the mass loss scheme; 
if dust grains with adequate extinction cross-sections condense, 
the momentum gained from photons which are absorbed or scat- 
tered by these particles may be sufficient to overcome the grav- 
itational pull of the star and accelerate the grains outwards. The 
collisions between dust particles and the surrounding gas cause 
a transfer of momentum that triggers a general outflow. 



1 



S. Bladh and S. Hofner: Exploring wind-driving dust species in cool luminous giants 



expected to be the main wind-drivers in C-type AGB stars. Self- 
consistent radiation-hydrodynamical models (RHD models) re- 
produce well observed wind properties and spectra of C-type 
AGB stars (e.g.|Le Bertre & Winters|[T998l|Winters et al.| |20lt)l 



Andersen et al. 2003 Gautschy-Loidl et al.| |2004| |Nowotny| 
etal.||2010||2011 |Sacuto et al.||201 1) , which demonstrates con- 
vincingly that the outflows are driven by radiative acceleration 
of amorphous carbon grains. 

For M-type AGB stars, on the other hand, where the dust 
chemistry is much more complex, it is still a matter of de- 
bate which grain species are responsible for driving the out- 
flows (see, e.g., the discussion in Hofner 2009). Observations 
of mid-infrared spectral features indicate that magnesium-iron 
silicates (olivine [Mg,Fe]2SiC<4 and pyroxene [Mg,Fe]SiC»3) are 
abundant in the circumstellar environment (e.g. |Molster et al.| 
|2010| l. However, the Mg/Fe ratio of the grain material is poorly 
constrained by such observations. The grain composition is im- 
portant because it strongly influences the absorption properties 
in the near-infrared, where the star emits most of its radiation, 
and consequently affects the resulting radiative acceleration. Fe- 
free silicates are extremely transparent in the near-infrared but 
the absorption cross-section increases dramatically with Fe con- 
tent. 

Chemical modeling based on phase-equilibrium between gas 
and dust favors silicate grains at the Mg-rich end of the spectrum, 
but the assumption of chemical equilibrium is not justified in 
the highly dynamical atmospheres of AGB stars. Kinetic mod- 
els of grain growth, on the other hand, point towards roughly 
equal amounts of Mg and Fe, reflecting the relative abundances 
of these elements (see, e.g., Gail 2010 for a detailed discus- 



sion on modeling of dust formation). The complex interplay of 
the stellar radiation field with the grain growth process, how- 
ever, has a decisive influence on the grain composition since the 
grain temperature at a given distance from the stellar surface is 
strongly affected by the Fe content. 

Using frequency-dependent wind models with a detailed 
treatment of dust formation, Woitke (2006) demonstrated that 
silicate grains have to be virtually Fe-free at distances corre- 
sponding to the wind acceleration zone. But the lack of Fe leads 
to low absorption cross-sections in the near-infrared region, re- 
sulting in insufficient radiative acceleration to produce outflows. 
Hofner (2008 ) suggested scattering as a possible solution: if con- 
ditions allow Fe-free silicate grains to grow into a size range 
of about 0.1 - 1/L/m, the scattering cross-section becomes dom- 
inant over absorption by several orders of magnitude, opening 
up the possibility of stellar winds driven by photon scattering. 
The models by Hofner (2008), which include a detailed descrip- 
tion for the growth of Mg2SiC>4 grains, show mass loss rates and 
wind velocities typical of AGB stars, as well as realistic photo- 
metric properties (cf. |Bladh et~aL 201 1 1. Strong observational 



support for the scenario of winds driven by photon scattering 
comes from recent work by |Norris et al.| ( |2012[ ), who detected 
dust particles of sizes ~ 0.3 /mi in the close circumstellar en- 
vironment of three M-type AGB stars, usingmulti-wavelength 
aperture-masking polarimetric interferometry^This answers the 
question whether grains can grow sufficiently large in the ex- 
tended atmosphere. 



1 The fraction of light scattered by dust is modeled assuming a geo- 
metrically thin shell of Fe-free silicates (forsterite) and a uniform grain 
size. The wavelength-dependence of the scattering cross-section allows 
them to derive the grain size from their multi-wavelengths measure- 
ments. 



Fe-free silicates are clearly a possible candidate for driving 
winds in M-type AGB stars but the question remains if other dust 
species can trigger winds or contribute to accelerating outflows. 
At present it is not clear whether winds are driven by a single 
type of grain material, a combination of several dust species or 
composite grains consisting of mixed materials (i.e. dirty grains 
or core mantel grains). In addition to silicates, observations show 
features attributed to other dust species such as spinel, corundum 
and simple oxides. 

Previous theoretical work on dust chemistry in AGB stars 
has focused on detailed descriptions of grain growth, investi- 
gating the time-dependent formation of homogenous and mixed 
grains, as well as networks of dust species competing for the 
same raw materials (e.g. Gail & Sedlmayr 1988 ] |Gauger et al. 



1990{|Gail & Sedrmayr||1999||Ferrarotti & Gail||2001[ Jeong 



etaL}|2003||Helling & Woitke , 2006; Woitke, 2006). Most exist- 
ing studies on dust-driven winds of M-type AGB stars (e.gpeong 
et al.| |2003| |Ferrarotti & Gail| |2006| l, however, are based on an 
crude treatment of the radiation field, resulting in unrealistic gas 
and grain temperatures and potentially misleading conclusions 
about the dust composition, as illustrated by the example of the 
Mg/Fe-ratio for silicate grains discussed above. 

The work presented here is part of an ongoing effort to in- 
vestigate how a strong radiation field affects dust condensation. 
To facilitate such an investigation we consciously scale back on 
the complexity of the dust properties and the grain growth pro- 
cess by using a parameterized formula for the dust opacity in 
our frequency-dependent RHD models. This formula captures 
the essential optical and chemical properties of the grain mate- 
rial when considering the interplay between the dust component 
and the radiation field. By constructing a grid of models where 
we systematically vary these properties, instead of conducting 
detailed studies for a few selected dust species, we get an overall 
picture of how the interaction between the radiation field and the 
dust component affects the dynamical properties of the model 
atmospheres. This knowledge makes it possible to exclude dust 
species which are not viable as wind-drivers and to identify pos- 
sible candidates for which a more detailed modeling is required. 

In the current paper we focus on the dynamics of the at- 
mospheres and winds. In a forthcoming paper we will give ad- 
ditional constraints based on comparisons of synthetic and ob- 
served spectra and photometric data. The paper is organized in 
the following way: In Sect. [2] we introduce simple analytical es- 
timates to pinpoint the criteria that have to be fulfilled for a dust 
species to be considered as a potential wind-driver. Sect. [3] in- 
cludes a description of the RHD wind model and the parame- 
terized dust opacity. In Sect.|4]we present the setup of the com- 
putational grid and how the input parameters were chosen. A 
presentation of the numerical results is given in Sect.|5]and we 
comment on specific dust species in Sect. [6] In Sect. |7]we pro- 
vide a summary of our conclusions. 



2. Dust driven winds: a qualitative picture 

Given the importance of dust in the mass loss mechanism of 
AGB stars it is essential to understand how the close stellar 
environment affects dust formation processes, and which grain 
properties are crucial to produce a mass outflow from the star. 
Some simple analytical constructs can help us gain qualitative 
insights into what properties of a grain material determine if a 
dust species is a possible wind-driver or not, paving the way for 
a physical interpretation of detailed RHD models. 
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2.1. Radiative acceleration of the wind 

Consider the dynamics of a fluid element, starting right after it 
has been accelerated outwards by a shock wave: the main forces 
acting on the fluid are the gravitational force and the radiative 
acceleration on the dust particles within the fluid element. The 
force due to the thermal pressure gradient in the atmosphere is 
small by comparison and we neglect it in this simple construct. 
It follows that the acceleration of the fluid element in the radial 
direction is mainly controlled by the relative value of the gravi- 
tational and radiative acceleration, which leads to the following 
equation of motion, 



du 
dt 



2grav + flrad = -Agrav (1 ~ D = ~g* I — j Q ~ 



(1) 



assuming a direct coupling between the gas and the dust com- 
ponents. Here u = dr/dt is the radial velocity of the fluid ele- 
ment at distance r from the stellar center, g* = GM„jR\ denotes 
the gravitational acceleration at the stellar radius R t and T is a 
dimensionless quantity, dependent on r, that measures the ratio 
between radiative and gravitational acceleration: 



flrad 



'grav 



4ncGM t ' 



(2) 



Included in the expression for F is the gravitational constant G, 
the speed of light c, the stellar mass and luminosity, M* and L t 
respectively, and the total flux-averaged dust opacity (k)h, 



<*>H 



So F * dA 



(3) 



where is the monochromatic flux at wavelength A and /c acc 
is the monochromatic grain opacity (see below). We define a 
critical value of the flux-averaged opacity when the gravitational 
and radiative acceleration are in balance, i.e. when r = 1. 



4ncGM t 



(4) 



Clearly the fluid element will accelerate outwards if F > 1 (i.e. if 
(k)h > K C iit) an d a closer examination of this quantity can help us 
pinpoint the grain properties necessary for driving a wind. The 
expression for F can be factorized into two parts; one part that 
solely depends on stellar parameters and fundamental physical 
constants, L*/47icGM*, and one part, (k)h, that is grain material 
dependent. Assuming for simplicity that all grains in the fluid 
element have equal radii a„, their collective opacity per mass of 
stellar matter can be expressed as 



(5) 



faccW, flgr) = -flgr2acc(/l,a g r)"gr, 

where « gr is the number density of grains and the efficiency Q- dQC 
is defined as the ratio between the radiative and the geometri- 
cal cross-section of an individual grain, which can be computed 
from optical data using Mie theory (e.g. Bohren & Huffman 
1 1983). Introducing the atomic weight of the monomer A mOB (the 
basic building block of the grain material), the bulk density of the 
grain material p gr , abundances of the limiting element £ij m and of 
helium su e , as well as the degree of condensation of the limiting 
element f c (the fraction of the limiting element condensed into 
grains), the expression for the dust opacity can be reformulated 
as (see Appendix |A| for details) 



^accO^s flgr) 



Q RCC (A, a gr ) em 



4 Pgr 



s(l +4e He ) 



fc 



(6) 



In this context we define the limiting element as the first element 
that will be completely consumed due to the relative element 
abundances in the atmosphere, adjusted for the number of atoms 
contributing to the monomer, s (stoichiometric coefficient). For 
example, in the case of Mg2SiC>4, the least abundant element 
would be Si (for a solar composition) but since 2 Mg atoms are 
used for building one monomer (s = 2), Mg becomes the limit- 
ing element^] The factor A mon /p gr is a material-dependent con- 
stant and the optical properties of the dust particles are captured 
in the efficiency Q^, which includes contributions from both 
absorption and scattering 



Gacc = Gabs + (1 - gsca)6sc 



(7) 



(gsca is the asymmetry factor describing deviations from 
isotropic scattering). In the small particle limit, where the grain 
radius is much smaller than the relevant wavelengths, the ab- 
sorption and scattering efficiencies behave like <2ab 
Gsca K flgr> according to Mie theory (e.g. 



oc q sr a nd 



Bohren & Huffman 



|1983[ ). In this limit absorption dominates over scattering, imply 
ing that gacc ~ Gabs and consequently that the efficiency per 
grain radius, 2acc/ fl gr ~ Gabs/ fl gr> becomes independent of grain 
size. We can therefore express the grain opacity during the early 
stages of grain formation, or during the whole process if parti- 
cles remain small, as a purely wavelength-dependent function, 



3A n 



4p g 



GaccW ' 



em 



i(l +4e He ) 



■/ c (27ra gr «,i) (8) 



where Q'. dcc = (2abs/flgr- Looking at each of the factors in this 
expression individually reveals what conditions need to be sat- 
isfied for a dust species to be a potential wind-driver. First, the 
efficiency per grain radius Q' acc has to be sufficiently large in the 
wavelength region of the stellar flux maximum, given that (k)h 
is calculated by taking the flux-mean of the grain opacity /c acc . 
Furthermore, the abundance £i; m is a limiting factor in the dust 
formation process. The last factor in the expression for the grain 
opacity, the degree of condensation f c , is a measure of how much 
of the limiting element has actually condensed into solid mate- 
rial. The degree of condensation will remain zero until the grains 
start to condense and grow. For this to happen, the gas has to be 
sufficiently cool and the grains have to be thermally stable. The 
grains are heated by the stellar radiation and consequently their 
temperature depends on both the optical properties of the grain 
material and the stellar flux distribution, as well as the distance 
from the star. The grain formation process is therefore strongly 
affected by the radiation field in the close stellar environment 
and the effect on grain temperature has to be taken into account 
when considering possible wind-driving dust species. 

2.2. Condensation and levitation distances 

In order to estimate if a specific grain material can start to con- 
dense in the density-enhanced layers of the atmosphere, we in- 
troduce the concepts of levitation distance and condensation dis- 
tance, and some easy-to-use approximations. In the following, 
the levitation distance Rf is defined as the distance to which 
pulsation-induced shock waves levitate the gas in the atmo- 
sphere, without radiative acceleration on dust. The condensa- 
tion distance R c is defined as the closest distance to the star 



2 Note that this is a simple way of estimating the maximum amount 
of condensable material for a specific dust species, and consequently, 
an upper limit for the opacity. This, however, does not necessarily mean 
that the addition of this element corresponds to the slowest rate (bottle- 
neck) in building up a monomer. 



3 



S. Bladh and S. Hofner: Exploring wind-driving dust species in cool luminous giants 

Table 1. Properties of a few selected dust species. Columns 2-7 list grain material properties, i.e. the atomic weight of the monomer 
A mon , the bulk density p gr , the limiting element and the adjusted abundance of said element s\\ m /s, the power law coefficient p and 
the condensation temperature T c . Column 8-10 include properties that depend on the chosen stellar parameters (see Sect. [4]), i.e. the 
condensation distance R c , the flux-averaged dust opacity for full condensation of the limiting element (/e)™ x and T. 



Material 


A 

^mon 


Pgr 


lim. 




P 


T c 


R c 


(cnr/g) 




r 




References & comments 




(«) 


(g/cm 3 ) 


element 






(K) 


(R,) 










Fe 


56 


7.87 


Fe 


3.24 • 10 -5 


2.4 


1050 


11.5 


1.5 




0.6 




1,4. 


A1 2 3 


102 


3.97 


Al 


1.48 • 10~ 6 


1.7 


1400 


3.6 


4 • 10~ 2 


2 


■ 10" 


-2 


1,6, 10. 


Ti0 2 


80 


4.26 


Ti 


9.77 • 10~ 8 


0.1 


1100 


3.4 


6 ■ 10~ 5 


2 


■ 10" 


-5 


2,3. 


Si0 2 


60 


2.20 


Si 


3.55 ■ 10~ 5 


0.2 


1050 


3.9 


6 • 10~ 2 


2 


• 10" 


-2 


1, 9. Interpolated optical data. 
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Refe rences. (l)prn1pO"To); (2)|Lattimer & Grossman|([T978l>; (3)|Zeidler et al.||20TT] >; (4) |Ordal et aLl([T988>; (5)|Rouleau & Martini ( [19911 ; (6) 
|Koike et al.|fl995| >; (7) pger et al.| l |2003) : (8) |Dorschner et al.|(T995> : (9) |Palik| < |T985> : (10) |Begemann et al.|fT997) ~ 

Notes. For the abundances of the limiting elements we adopted the values from Grevesse & Anders ( 19891, except for C, N and O where we took 
the data from Grevesse & Sauval ( 1994). For amC the abundance was calculated using the formula e c = (C/O - l)so, with a C/O-ratio of 1.25. 
The first reference in column 1 1 corresponds to condensation temperature (from chemical equilibrium calculations) and the second to bulk density 
and optical data in the NIR. The values given for Si0 2 in the NIR are derived from an interpolation, using data from Palik ( 1985 1, and probably 
overestimating the cross-sections. SPL is short for small particle limit. 
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Fig. 1. Curves of constant condensation distance (R c /R* = 
2,4, 10) as a function of the power law coefficient p and the 



condensation temperature T c , using Eq. ( 10 1, with T* = 2800 K 



(solid lines) and T t = 2500 K (dashed lines). The filled red cir- 
cles mark where the selected set of dust species are situated in 
the /?/r c -plane. The dust species in parenthesis have uncertain or 
interpolated optical data in the near-infrared. For T c and p values 
for the individual dust species, see Tab. [T] 



were grains of a specific type can exist, i.e. are thermally sta- 
ble. For a grain material to be considered a possible wind-driver, 
the two stages of the mass-loss scheme have to connect. The 
second stage, i.e. the radiative acceleration, can only be initiated 
if levitation by shock waves lifts gas beyond the condensation 
distance. 

A simple argument using the complete conversion of 
pulsation-induced kinetic energy into potential energy, ignoring 



heat loss and pressure effects, can give us an estimate of how 
high shock waves can lift gas. If we assume the gas has an ini- 
tial velocity uq at distance Ro from the center of the star we can 
derive the following expression for the levitation distance 



Re _ Ro 
R~*~ R~ t 



Ro / uq 



R* \M, 



(9) 



where the escape velocity at the stellar surface is given by 
«esc = (2M i G/R„) 1 ^ 2 (see Appendix O. Given stellar parameters 
typical of an AGB star, i.e. M, = 1 M Q , L„ = 5000 L Q and T eff = 
2800 K, the escape velocity is about 36 km/s. Radial veloci- 
ties derived from observations of the second overtone CO line 
(Av = 3) are of the order of 10-15 km/s. According to dynam- 
ical models, these lines are formed in the region from the stel- 
lar surface out to about 1.5 stellar radii (Nowotny et al. 2010). 
Assuming an initial velocity of uq = 15 km/s and a distance 
Ro = 1.5 R m , Eq. ([9]) gives a levitation distance of about 2R t . 
This corresponds approximately to where interferometric mea- 
surements place the inner edges of dust shells around AGB stars 
(e.g. 



et al. 



Wittkowski et al. 
|2012) . 



2007 Karovicova et al. 201 1 Norris 



Assuming a Planckian radiation field, geometrically diluted 
with distance from the star (i.e. an optically thin atmosphere 
with no molecular features), a power law for the dust absorp- 
tion coefficient /c a b s in the relevant wavelength range, and that 
the grain temperature is determined by the condition of radiative 
equilibrium, the condensation distance can be expressed as (see 
Appendix |B) 



R t 



1 (To 

2 \T t 



where oc A p , 



(10) 



where T* is the effective temperature of the star. In Fig. [T] we 
have plotted curves of constant condensation distances as a func- 
tion of p and T c , using this formula. The chosen values of effec- 
tive temperature, T* = 2800 K and T* = 2500 K, correspond 
to the detailed models used in this paper and by Woitke (2006 ), 
respectively. Due to the wavelength dependence of the absorp- 
tion coefficient (i.e. the value of p) different condensates will 
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react differently to the close stellar environment. A grain mate- 
rial with a positive value of p tends to heat up when interacting 
with the stellar radiation field, being more efficient at absorbing 
than emitting radiation, and thereby pushing the condensation 
distance further out. The opposite is true for a grain material 
with a negative value of p. To summarize, it is the combina- 
tion of condensation temperature and the slope of the absorption 
coefficient that determines how close to a star a grain type can 
survive. For example, Fe-bearing silicates have a condensation 
temperature of about HOOK and p m 2.3, which according to 
this simple formula results in a condensation distance of around 
10R* for T* = 2800 K, far outside the reach of the levitated gas. 
The corresponding value for Fe-free silicates are T c as HOOK 
and p as -0.9, resulting in a condensation distance of 2R* for 
T* = 2800 K, which is within reach of the levitated gas. 

As mentioned in the introduction, phase-equilibrium cannot 
properly describe the dust formation in AGB-stars. According 
to detailed models, the formation process demands cooling well 
below the condensation temperature to initiate grain growth (see, 
e.g., Gail 20 1 0| > . Therefore, assuming thermal stability will lead 
to an underestimation of the condensation distance, and we treat 
it here as a lower limit. A 10% decrease of the condensation 
temperature in Eq. ( 10 1 leads to an increase in condensation dis- 
depending on p, where a positive value 
Despite the simplifications 



tance between 17-41% 
of p results in larger changes in R, 
involved in Eq. ( 10 1, however, it turns out that the resulting con- 
densation distances for amorphous carbon and Fe-free silicates 
compare well to RHD models which combine a detailed descrip- 
tion of grain growth with a frequency-dependent treatment of ra- 
diative transfer (e.g. Hofner et al. 2003 Woitke 2006; Hofner 
2008 ). These models place the condensation zone of such grains 
typically in the range 2-3/?*. 

2.3. Criteria for wind-drivers 

To summarize, we have the following properties that determine 
if a specific grain material (pure or mixed) can be excluded as a 
wind-driver or be regarded as a potential candidate, where fur- 
ther investigation is needed: 

(a) The distance from the star where the grains are thermally 
stable, determined by the condensation temperature T c and 
the near-infrared slope of the absorption coefficient p. 

(b) The abundance of the limiting element of the grain material. 

(c) The absorption and/or scattering efficiency of the grain ma- 
terial in the wavelength region where most of the stellar ra- 
diation is emitted. 

Further criteria, such as time-scales associated with growth rates 
that determine if a dust species can grow fast enough, have to be 
investigated with detailed non-equilibrium models of the atmo- 
spheric chemistry (for a more in depth discussion on this topic 
see, e.g., |Gail| [2010 or the references given in the introduction). 

Estimates for specific grain materials can be obtained in the 
following way: if the particles can form sufficiently close to the 
stellar surface, as can be checked by Eq. Q and Eq. ( |T0] > (levita- 
tion distance and condensation distance, respectively), the com- 
bined effect of the critical abundance s\i m and the efficiency Q- dcc 
can be investigated by setting f c = 1 in the expression for the 
dust opacity /e acc (Eq. ()8]l) and assuming a suitable flux distri- 
bution. If the resulting flux-averaged dust opacity (k)h is larger 
than the critical opacity A: cr i t the grain material is a possible can- 
didate for triggering outflows. Examples of (k) h for different 
dust species, using the flux distribution plotted in the bottom 
panel of Fig. [2] (black curve), are given in Tab. [T] If one or more 
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Fig. 2. The top panel shows the efficiency per grain radius 
Gacc/flgr, in the small particle limit, as a function of wavelength 
for a selection of dust species. The dashed lines show the power 
law fit according to Eq. ( (20) ) and the shaded area indicates the 
wavelength region for which the optical data is fitted. The nor- 
malized flux of the hydrostatic initial model (black) and the 
Planck function for T e g = 2800 K (red) are plotted in the lower 
panel. Due to absorption by TiO molecules in the visual the flux 
of the initial model is very non-Planckian in this wavelength re- 
gion. 

of the above criteria is not fulfilled for a specific grain material 
then we can exclude it as a wind-driver since the analytical ap- 
proximations assume the most favorable limit. 



3. Detailed RHD models 

The analytical estimates discussed above presents a very sim- 
plified qualitative picture of the physics in the atmosphere of 
AGB stars; the dynamical effects of the gas pressure are ig- 
nored, the radiation field is described by a geometrically di- 
luted Planck function and F is considered as a simple function of 
the flux-averaged grain opacity. In reality, the molecular opaci- 
ties strongly affect the radiation field, pulsation-induced shock 
waves influence the structure of the atmosphere and F is a time- 
varying function depending on grain growth. To be able to take 
into account the complex interplay between gas dynamics and 
the radiation field we use state-of-the-art RHD models (Hofner 
|et al.[ 2003 ; Hofner, 2008 and references therein) that have been 



tested against different types of observations (e.g. Gautschy 
ILoidl et all [2004) |Nowofny et al.| [2010) |20H~) |Sacuto et al 



|201 l)|Bladh et al.f|201 l| i, and can provide us with realistic tern 
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perature and density structures. In these numerical models we 
include a simple parameterized description of the dust opacity, 
designed to investigate what dust species can exist close to the 
stellar surface and potentially drive outflows. 

3.1. Gas dynamics and radiation field 

The model atmosphere covers a spherical shell with an in- 
ner boundary situated just below the photosphere and an outer 
boundary in accordance with the dynamical properties of the 
model. In models that develop winds the outer boundary is fixed 
at the point where the flow velocity has reached its terminal 
value, allowing outflow, and in models without winds the outer 
boundary follows the periodic motions of the upper atmospheric 
layers. The variable structure of the atmosphere is described by 
the equations of hydrodynamics (equation of continuity, equa- 
tion of motion and the energy equation) and the pulsations are 
simulated by temporal variation of physical quantities at the in- 
ner boundary. The opacities of molecules and dust forming in the 
outer cool layers of the atmosphere strongly affect the radiation 
field and in order to achieve realistic density and temperature 
structures the models include a frequency-dependent treatment 
of the radiative transfer (see Ho fner et al.[ |2003[ for more de- 
tails). 

The conservation of mass, momentum and energy is de- 
scribed by the following equations for the gas component: 



simultaneously with Eq. (TT)-([T3|. Note that the terms contain- 
ing the dust opacity cancel out in Eq. ( p~5] > since we have assumed 



that radiative equilibrium holds for the dust grains (see Eq. (|14|) 
The remaining unknown quantities required for closing the 
system of conservation laws (i.e. the frequency-averaged opaci- 
ties, Aj, kh and and the Eddington factor _f e dd = K/J) are de- 
termined in a separate step: by solving the frequency-dependent 
radiative transfer equation for the current density-temperature 
structure after each hydrodynamical time-step we obtain J v , H v 
and K v , which allows us to compute the averaged opacities 



i r°° 



Xydv where X 



-r 

Jo 



Xydv. 



(17) 



dt 
a 

at 
a 

at 



(p) + V • (pu) = 



(11) 



and the Eddington factor that are used in Eq. ( 12 H 16 1. In the 
current models we use 319 wavelength points, in contrast to 51 
points in Hofn e~t al.| p003), for improved representation of the 
opacities and the radiation field. Assuming LTE, we approximate 
the source functions with Planck functions, i.e. S s = B V (T S ) 
and Sd — ByiJi), and the corresponding averaged opacities are 
Planck means. 

It should be mentioned here that we neglect scattering when 
solving the frequency-dependent radiative transfer equation be- 
tween hydro steps. Scattering on dust particles can have a sig- 
nificant effect on the momentum gain of individual grains and 
thus the wind dynamics. However, scattered photons will not 
noticeably change the overall radiation field, especially as the 
grains grow beyond the small particle limit and scattering be- 
comes more forward oriented. Therefore the effect of scattering 



/ \ . / \ {jvn r 47rp / g j\ D .... on the frequency-dependent intensity and its moments J y , H v , 

(pu) + v • (puu) = -\Po *-p + [Kt+K U )H (12) . . J . . v ' v ' 

° r c v ' and K.. which arp rpnnirpn rn rnmnnte \ne avprasTP nnanrips (spp 

(pe) + V ■ (peu) = -P g V -u+Anp (k]J - k^S g ) (13) 



where u is the radial velocity, p is the gas density, P g (p, e) is the 
thermal gas pressure, m r is the integrated mass within a spher- 
ical shell with radius r and e is the specific internal energy of 
the gas. We assume direct coupling between the motion of the 
gas and the dust, i.e. the momentum gained by the dust from 
the radiation field is directly transferred to the gas. Both the ab- 
sorption and the scattering efficiency of the grain material may 
contribute to the overall momentum gain (see Eq. (j7]i) and both 



are included in In Eq. (|12|) and d 1 3 b we use a simplified no 



tation for the opacities compared toThe rest of this paper. The 
superscripts g and d correspond to the gas and dust components 
respectively, and the subscripts indicate averages over different 
moments of the radiative intensity (see below). This means that 
4 corresponds to (k)h as defined by Eq. (Hi, (|6ji and d7|. 

The energy budget of the dust component, and therefore the 
grain temperature, is determined by the condition of radiative 
equilibrium 



\ibs,S 



S(T£ = 0—>T d 



abs,J 



V"abs,S/ 



1/4 



(14) 



In this formula T r = y/Jir/cr denotes the radiation tempera- 
ture and /e a bs the true absorption part of the dust opacity. The 
frequency-integrated moments of the intensity, J and H, are de- 
termined by solving the zeroth and first moment equation of the 
radiative transfer equation 



and ICy, which are required to compute the average opacities (see 
Eq. ( 17 1), will be smallj^The exclusion of scattering in the radia- 
tive transfer leads to a significant reduction of the computational 
effort and allows us to use a larger number of frequency points, 
which is important for obtaining realistic temperature structures. 

3.2. Parameterized dust description 

To systematically study the effect of a close stellar environment 
on different grain materials we construct a simple parameterized 
formula describing the dust opacity. This approach allows us to 
simulate aspects such as composite grain materials and effects of 
absorptions versus scattering cross-sections in a simple way, in 
addition to circumventing the problem of missing or incomplete 
optical data for certain dust species in the near infrared wave- 
length region (cf. Zeidl eret al.||201 !"] ). Since we are interested 
in dust formation in the vicinity of a strong radiation source we 
focus on the factors in the expression for the grain opacity for 
which the radiation field will be important. One such factor is 
how the efficiency per grain radius, Qacc/ctgr, varies with wave- 
length, which in combination with the radiation field determines 
the grain temperature. Another important factor is the condensa- 
tion temperature, a chemical property that indicates below which 
temperature grains of a certain material are thermally stable. We 
therefore construct a parameterized dust description that allows 
for both different optical wavelength dependences and conden- 
sation temperatures, i.e., 



c (A)=k(A)-f c (r,t,T c ). 



(18) 



V -H+p(i<$J-K*S e ) = 
3K — J 

VK + — — +p(4+4)// = o 



(15) 
(16) 



3 Note that scattering is included in the frequency-dependent opaci- 
ties k v in Eq. as appropriate for the respective means, i.e. k v cor- 
responds to k. xc when calculating the flux mean (k)^ that enters the 
equation of motion Eq. 1 12 1. 
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Fig. 3. The degree of condensation as a function of grain tem- 
perature, using the parameterized dust description outlined in 
Sect|3.2| (solid black line), and for a model with a detailed de- 



scription of Mg2Si04 grains (Hofner 2008 model A; dashed 



lines). The different dashed curves represent a selection of pul- 
sation phases ip from the detailed model, where the degree of 
condensation has been normalized to vary between and 1 to 
facilitate comparing the shape of the curves. The grey zone indi- 
cates how the parameter Ar sets the width of the dust formation 
zone. 



Our description is inspired by a formula used by Bowen ( 1988 ) 
in his dust driven wind models, but is here generalized to allow 
for wavelength-dependent optical properties. The degree of con- 
densation, / c (r, t, T c ), is designed to increase monotonically with 
falling grain temperature, approaching a value of 1 as the grain 
temperature drops well below the condensation temperature (see 
Fig# 



fc(r,t,T c ) = 



1 



1 + e (Ta(r,t)-T c )IKT 



(19) 



The grain temperature Ta is determined by condition of radia- 
tive equilibrium (Eq. ( 14 1) and it is a function of both time and 
distance from the star, due to the varying radiation field. The pa- 
rameter T c sets where f c = 0.5 and Ar regulates the width of the 
dust formation zone. The value of AT is chosen in accordance 
with models that include detailed dust formation (see Sect. [4j. 
The wavelength-dependent part of the dust opacity k is modeled 
as 



k(A) = kq 



(20) 



where p is obtained by fitting a power-law function to 2abs/flgi 
data in the wavelength region where most of the stellar flux is 
emitted (see Fig. [2] and Sect. |4jl. Regions of low flux will not 
contribute significantly to the energy balance and the radiative 
acceleration. In this expression kq is a scaling factor such that 
k(Ao) = kq. In general the dust opacity would be a function of 
both wavelength and grain radius. However, since we are inter- 
ested in the onset of grain formation and in particular the lo- 
cation of the dust formation zone, the adopted parameterization 
mimics the small particle limit where QaccA*gr ~ 6abs/flgr is not 
a function of grain size (small particles have to exist before they 
can grow beyond the small particle limit). 

To distinguish between effects of scattering and absorption 
of photons on dust grains we further introduce a quantity / a b s 



which sets the percentage of the dust opacity /c a , 
considered as true absorption 



*abs(/0 = /abs • k(X) ■ f c (r, t) where / at 



^abs 



that is to be 



(21) 



The dust opacity K acc , which includes contributions from both the 
absorption and scattering cross-sections, is used to calculate the 
radiative acceleration in the equation of motion, and the true ab- 
sorption part /c a bs is used to determine the grain temperature (see 
Sect. - 



3.1 



This allows us to separate the dynamical and thermal 
effects of the dust opacity and the parameter / a b s can be adjusted 
to explore the effects of varying degrees of grain transparency. 
This parameter can also be used to simulate additional forces 
that affect the dynamics of the gas without changing the energy 
distribution of the radiation field. 



4. Model parameters 

The stellar parameters chosen for the dynamical models in this 
study are typical for an M-type AGB star, with a stellar mass 
and luminosity of 1 M Q and 5000 L Q , respectively, an effective 
temperature of 2800 K and solar abundances. The velocity am- 
plitude used to simulate pulsations at the inner boundary was 
set to 4 km/s, which has proven to be a reasonable value when 
comparing with observations of molecular lines forming in the 
inner atmospheric region ( |Nowotny et al.||2~010| l. The same stel- 
lar parameters and piston velocity were used by |Hofne"rj ( 2008) 
in models with a detailed grain growth description for Mg2SiC>4 
particles, demonstrating that such grains can drive winds if they 
grow to sizes where scattering dominates the dust opacity. 

Our formula for the dust opacity contains the optical param- 
eters p and kq, as well as the condensation temperature T c and 
the width of the condensation zone Ar. For the purpose of this 
study p and T c will be considered as free parameters within a 
reasonable range. We focus on how the grain temperature, and 
consequently the condensation distance, is affected by the stellar 
radiation field (criterion (a) in Sect. |2.3) , disregarding the possi- 
bility of too low abundances or cross-sections for specific grain 
species (criteria (b) and (c)) in our grid of RHD models. We 
therefore choose a value for a:o such that Y > 1 when the grains 
have fully condensed (f c = 1), making sure that the dust opacity 
will be high enough to trigger an outflow if grains condense suf- 
ficiently close to the star. For our chosen set of stellar parameters 
we find /c cr ; t = 2.6 from Eq. (|4]). In the parameterized dust opac- 
ity we set kq - 3.0 when p — and then adjust for other values 
of p so that (k)h remains fixed. To provide good estimates for 
kq we flux-average over the spectrum of the hydrostatic initial 
model instead of a Planckian radiation field, since absorption by 
TiO molecules in the visual makes the flux distribution very non- 
Planckian. This results in a flux-averaged opacity (k) h m 3.0 and 
a corresponding value Fs 1.1, for all values of p in the dynam- 
ical models. 

The parameter Ar that adjusts the width of the dust forma- 
tion zone in the parameterized dust description is set to 100 K. 
This value was chosen by comparing with the typical width of 
the dust formation zone in a model with a detailed dust descrip- 
tion and the same stellar parameters (see Fig. [3}. The exact value 
of Ar does not have a significant effect on the resulting wind 
properties. 

Finally, with the stellar and dust parameters chosen, we set 
up a grid to investigate the dynamical effects of different com- 
binations of chemical and optical properties in potential wind- 
driving dust species by varying the input parameters p and T c in 
the dust opacity. To cover the most probable combinations we 
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Fig. 4. The panels show wind velocities, mass loss rates, radiative acceleration and density for a grid of dynamical models with 
stellar parameters M, = 1 M Q , L* = 5000 L G and T e g = 2800 K, a piston velocity of 4 km/s and / a b s = 1 .0. The grid covers a range 
of chemical and optical dust properties, represented by the variables T c and p, and includes in total 165 dynamical models. The 
sawtooth-like shape of the boundary between models with and without a wind is a consequence of the discrete number of models in 
the grid. The radiative acceleration and density is measured at the distance from the star when f c = 0.5 (where 7^ = T c ), averaged 
over pulsation cycle, and the velocity and mass loss rate are measured at the outer boundary of the model. The over-plotted contours 
are curves of constant condensation distance (R c = 2,4, 10) using Eq. ( fTO] ). 



let p vary from -1 to 2.5 and T c from 700K up to 1700 K, in 
increments of 0.25 and 100K respectively, including a total of 
165 sampling points (individual dynamical models). In addition, 
we also vary the degree to which the dust opacity is considered 
true absorption by setting varying the parameter / a b s to 1.0 and 
0.5 respectively (100% and 50% true absorption). 



5. Results 

The detailed dynamical models provide us with the radial struc- 
ture of the atmosphere and wind, from the photosphere to where 
the outflow reaches its terminal velocity (typically 20 - 30 R*), 
as a function of pulsation phase. This includes velocities, 
temperature- and density-structures and degree of condensation 
as a function of distance for the dust component. The mass loss 
rate and wind velocity at the outer boundary are recorded at each 



time-step and the values calculated for individual models are av- 
eraged over about 500 periods. 

5.1. Constraints on potential wind-drivers 

The grid of dynamical models shows very distinct trends as to 
what combination of optical and chemical dust properties, i.e. T c 
and p, are necessary for driving a wind in M-type AGB stars (see 
Fig.Q. Dust material with high condensation temperatures and 
a NIR absorption coefficient that increases, or decreases slowly, 
with wavelength will condense close enough to the stellar sur- 
face to be able to trigger an outflow. As an example, dust mate- 
rials with a condensation temperature around 1700 K and a neg- 
ative slope as steep as p * 1.5 still produce winds. However, 
for condensation temperatures below 1400 K we need p < to 
trigger an outflow. 
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Fig. 5. Dust-, gas- and radiation-temperatures as a function of 
distance for three different models where the condensation tem- 
perature and phase are kept fixed (r c = 1500 K and tp = 0.25) 
but allowing for different values of p. The dashed lines show the 
temperature profiles of the hydrostatic starting model for the gas 
component (black) and the radiation field (blue). 



The shape of the boundary between the wind and no-wind 
regime resembles the curves of constant condensation distance 
according to the simple estimate in Eq. ( 10 1 (see the over-plotted 
contours in Fig. [4]). This demonstrates the crucial role of the 
condensation distance in connecting the two stages of the mass 
loss scheme. While there is a qualitative agreement between the 
trends in the detailed models and the simple estimates it should 
be noted that the latter will not give quantitatively correct results. 
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Fig. 6. Dust-, gas- and radiation-temperatures as a function of 
distance for a model with a time-dependent description of the 

The 



2008). 



dust formation, using Mg2Si04 grains (Hofner 
dashed lines show the temperature profiles of the hydrostatic 
starting model for the gas component (black) and the radiation 
field (blue) at <p = 0.25. 



This can be seen in Pig. [9j showing the actual condensation dis- 
tance in the detailed RHD modelsQ 

The exact location of the boundary between the wind and 
no-wind zone will depend on stellar parameters, pulsation am- 
plitudes and certain physical assumptions made in the model. 
In particular, the simple parameterized dust opacity will tend to 
over-estimate the efficiency of dust formation due to the instant 
coupling between the degree of condensation and the grain tem- 
perature. That means that the dust opacity will take effect as soon 
as the grain temperature drops below the condensation temper- 
ature and the time-scales for grain growth are ignored. To com- 
pensate for this effect, we only include models that show well- 
developed winds, eliminating a few cases where a more sophis- 
ticated treatment of the dust formation probably would fail to 
produce stable outflows. Regarding stellar parameters, we focus 
the discussion in this paper on one set of values which are typical 
for M-type AGB stars, due to the computational efforts involved 
in creating a grid of this size. As far as the condensation distance 
is concerned, the most important stellar parameter is the effec- 
tive temperature. We can estimate the effect of a change in effec- 
tive temperature by using Eq. ( 10 1. As can be seen in Fig.[T] the 
curves of constant condensation distance shift by about 100K to 
left when we change T* = 2800 K to T t = 2500 K. The bound- 
ary between the wind and no-wind zones is therefore probably 
not as sharp as depicted in Fig. [4] and [8] but more like a smooth 
transition zone. 



5.2. Trends in grain temperature 

The shape of the boundary between models with and without 
winds in the p/r c -plane (see Fig. [4]) is a consequence of the 
strong dependence of the grain temperature on p. A negative 
value of p corresponds to a positive slope of the efficiency per 
grain radius Q acc /a gr (see Fig.|2ji in the wavelength region where 
most of the stellar flux is concentrated and materials with such 
optical properties will emit radiation more easily than they ab- 
sorb. Therefore, in radiative equilibrium, the grain temperature 



4 Note in this context that the condensation distance is not a param- 
eter in the detailed RHD models, but a result of the complex interplay 
between the dust properties and the radiation field. 
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Fig. 7. Trends in velocity and mass loss rate as a function of T c for p varying betwen [-1, 1]. In the left two panels / a b s = 1.0 and in 
the right two panels / a b s = 0.5. 



will be lower than the radiation temperature, allowing the grains 
to form closer to the stellar surface. The opposite will happen 
for positive values of p, where the slope of Q acc /a gT is negative, 
and dust particles instead absorb radiation more easily than they 
emit, pushing the condensation distance further away from the 
stellar surface. This trend in grain temperature for positive ver- 
sus negative values of p is clearly visible in Fig. [5] where the 
temperature profiles for the gas component, the dust component 
and the radiation field are plotted as a function of distance for 
a fixed condensation temperature and phase, but with different 
values of p. In the top panel p = 1 and when dust starts to form 
the grain temperature becomes warmer than the radiation tem- 
perature. In the middle panel p is set to zero and consequently 
the radiation temperature and dust temperature are equal. In the 
bottom panel p = — 1 and we see a cooler dust temperature com- 
pared to the radiation temperature as soon as dust starts to form. 

For comparison we also plot the temperature structure for a 
model which includes a detailed description of Mg2SiC>4 grains 
regarding both their formation and optical properties (see Fig. [6}. 
Using the power law fit defined in Sect. 3.2 we obtain p as -0.9 
for this material. However, as can be seen in Fig 3, a power law 
is not a good representation of the data outside the marked grey 



area, underestimating the contribution to cooling at longer wave- 
lengths. As a consequence, the grain temperature in the detailed 
model is even lower than for the corresponding parameterized 
case (p = — 1, bottom panel of Fig. [5]). 



5.3. Trends in mass loss and wind velocity 

Since we chose kq to produce a flux-averaged dust opacity large 
enough to trigger a wind when grains have fully condensed 
(F > 1 for f c = 1), dynamical properties like mass loss rates 
and wind velocities will not be representative of specific dust 
species. Despite this, there is still interesting information to de- 
rive from the overall trends in mass loss rate and wind velocity 
in this parameter study. For the mass loss rates, as can be seen 
in the bottom panels of Fig. [7] there is a strong correlation be- 
tween increasing condensation temperature and increasing mass 
loss rates. There is also a weaker correlation between decreas- 
ing value of the power law coefficient p and increasing mass 
loss rates. Both increasing the condensation temperature and de- 
creasing the value of the power law coefficient pushes the con- 
densation distance closer to the stellar surface, where the density 
p is higher, which is what causes the higher mass loss rates since 
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M oc p. The correlation between mass loss rate and density at the 
condensation distance is evident in the bottom panels of Fig. [4] 

The radiative acceleration, and consequently the velocity, of 
the out-flowing gas depends on the flux-integrated dust opacity 
(k)h. This quantity is affected by both the spectral energy distri- 
bution of the stellar flux and the power law representation of the 
dust opacity (see Eq. (|3]l). As mentioned above, increasing the 
condensation temperature and decreasing the value of the power 
law coefficient pushes the condensation distance closer to the 
stellar surface and results in higher mass loss rates. As dust starts 
to form and the optical depth of the wind increases, the spectral 
distribution of the radiation field will shift towards longer wave- 
lengths compared to the photospheric energy distribution. The 
redding of the radiation field will be more pronounced the higher 
the mass loss rate is and the more opaque the dust particles are. 

The radiative acceleration of dust material with different 
power law coefficients will be affected differently by the redden- 
ing of the radiation field: a negative value of p (positive slope 
of 2/flg r ) will result in a higher radiative acceleration compared 
to a positive value of p (negative slope of Q/a gr ), due to the in- 
creased radiation in the far infrared. This gives rise an increasing 
spread in velocity for different values of p with increasing con- 
densation temperature, which can be seen in the top left panel in 
Fig. 7 This trend is less pronounced (see the top right panel in 
Fig. 7 1 when we assume more transparent grains by decreasing 
the value of / a b s without decreasing the radiative acceleration, 
since the reddening of the radiation field is reduced when the 
dust absorption is lowered. 



5.4. Effects of varying degree of grain transparency 

In addition to exploring the dynamical effects of different opti- 
cal wavelength dependences and condensation temperatures we 
also investigate the effects of grains that show various degrees of 
transparency. Both scattering and true absorption may contribute 
to the radiative acceleration of dust particles (see Sect |2. 1 1 > and 
scattering might even be the dominant process for momentum 
transfer, as demonstrated for Mg2SiC>4 grains by Hofner (2008 ). 
As in the case of the different values of p, we keep the radia- 
tive acceleration fixed but vary the fraction of K acc that is consid- 
ered as true absorption. For the values of / a b s that we have tested 
we obtain a similar pattern in the parameter space concerning 
the wind/no-wind zone; combinations of high positive ^-values 
and low condensation temperatures produce no outflows (see the 
right panels of Fig. |4] and [8]>. The location of the boundary be- 
tween the wind and the no- wind regime is, however, slightly 
shifted. This is a consequence of how the dust opacity affects 
the flux distribution. A higher degree of true absorption by the 
dust grains results in a stronger reddening and thermalization of 
the radiation field. The larger optical depth of the circumstel- 
lar envelope leads to a more Planckian radiation field. This may 
explain why the shape of the boundary in Fig. |4]follows the con- 
tours of constant R c derived from Eq. ( 10 1, which is based on a 



Condensation distance - [R« 



Planckian flux distribution, more closely. 

A detailed comparison between the models with / a b s = 0.5 
and / a b s = 1.0 shows somewhat different dynamical structures 
in the atmospheres that produce winds. The mass loss rate is 
strongly correlated with the density in the dust formation zone in 
both cases (see the bottom panels of Fig.|4]and[8]), but the radial 
position of the zone itself varies less for models with / a b s = 0.5 
(see Tab. [2). The smaller variation in velocity in the grid with 
/ a bs set to 0.5 is expected, since the circumstellar reddening of 
the radiation field is less pronounced when dust particles absorb 
less of the stellar flux. The variation in mass loss rates for the two 
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Fig. 9. Distance from the stellar surface where the degree of con- 
densation has reached f c = 0.5 (where T& = T c ), averaged over 
pulsation cycle. In the top panel / a b s = 1.0 and in the bottom 
panel / abs = 0.5. 

Table 2. Ranges of dynamical and atmospheric properties in 
models producing outflows in the two RHD grids (/ a b s = 0.5 
and/ abs = 1.0). 







/abs 


= 0.5 


/abs 


= 1.0 


u 


[km/s] 


9.5 


- 16.0 


5.9 


-20.5 


M 


[Afo/yr] 


i • i(r 7 


- 3 • 10- 6 


7 • 10- 8 


- 3 • 10- 6 


R c 




2.3 


-2.5 


2.2 


-2.8 


P(T C ) 


[g/cm 2 ] 


4- i(r 15 


-4- 10" 14 


6 • 10~ 16 


- 5 ■ 10" 14 



sets of models is a consequence of the variation in the density at 
the condensation distance. 



6. Discussion: specific grain materials 

To demonstrate how the result presented in the previous sections 
can be used to evaluate potential wind-drivers, let us take a closer 
look at a few specific grain materials relevant for M-type AGB 
stars. Typical values of p can be derived from optical data of ac- 
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Fig. 8. Same as Fig.^ but for / abs = 0.5. 
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tual grain materials, as shown in Fig. |]and Tab. _ 
metallic iron, Fe-bearing silicates and amorphous carbon can all 
be fitted well by a power-law function, even quite a bit outside 
the wavelength range where the star radiates most of its flux. A 
power-law fit for the Fe-free silicate, however, is less appropriate 
outside the considered wavelength range and while this approach 
captures the essential dynamical effects of the dust material, care 
has to be taken when considering photometric fluxes outside the 
fitted region (Bladh et al. 2012 in prep.). 

The presence of metallic iron in the circumstellar environ- 
ment of AGB stars has been suggested based on observations 
(e.g|McDonald ~et~aL"||2010| l. This material is very opaque in the 
near-infrared wavelength region (see Pig.[2j. However, assuming 
that all available iron has condensed into dust and averaging over 
the photospheric flux distribution results in a total opacity which 
is lower than the critical opacity needed to drive a wind. More 
importantly, due to the combination of p as 2.4 and T c m 1050 K 
for this material, iron grains cannot form at distances that can be 



5 The wavelength region where the power law is fitted, indicated by 
the grey area in Fig. [2] corresponds to where the condition E > |0.5£ max | 



is satisfied. In this expression E = AB A (T,,A), where E lr 
and A max is given by Wien's law. 



E(Amax) 



reached by the shock-levitated gas and can therefore not trigger 
an outflow (see Fig.[T]and|4]l. 

A better suggestion for possible wind-drivers are Fe-bearing 
silicates (MgFeSiCU). They have an absorption efficiency com- 
parable to metallic iron in the near-infrared but are not as com- 
pact (lower density). This results in a flux-averaged dust opac- 
ity above the critical opacity for our chosen stellar parameters, 
when assuming full condensation. But as can be seen from Fig. 
[T] and |4j the optical and chemical properties of Fe-bearing sil- 
icates (p as 2.3 and T c ai HOOK) prevent them from forming 
sufficiently close to the stellar surface. 

In contrast, Fe-free silicates such as Mg2SiC>4 or MgSiC>3 
can form in close proximity to the stellar surface but are very 
transparent in the wavelength range around 1 /jm, leading to a 
total opacity well below the critical opacity as long at the parti- 
cles are small. But if the particles grow to grain sizes comparable 
to the wavelength of the flux maximum, the contribution to the 
efficiency Q acc from the scattering cross-section is substantial, 
and for the case of Mg2SiC>4 grains, scattering has been demon- 
strated to be sufficient to drive a wind (seepofner, 2008] Fig. 1). 
In the framework of our parameterized RHD models, this ma- 
terial has p as -0.9 and T c ai HOOK, which puts it right on 
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the border of the wind/no-wind zone in Fig. |4] However, as dis- 
cussed in Sect. 5.2 the power law fit tends to overestimate the 
grain temperature compared to the real optical data, resulting in 
a larger condensation distance than in a detailed model. 

The transparency of T1O2 together with the low abundance of 
the limiting element rules this material out as a potential wind- 
driver, even if it can form close to the stellar surface. Lastly, it 
is difficult to say anything definite about dust species like SiC>2 
or AI2O3 at present due to uncertainties in measured optical data 
and incomplete wavelength coverage. 

Optical properties are predominantly measured in wave- 
length regions where there are distinct spectral features and lab- 
oratory studies often leave out regions where opacities are ex- 



pected to be comparatively low (see discussion in Zeidler et al. 
201 1| . Characteristic features are important for identifying indi 
vidual grain materials but may not be essential for the overall 
energy balance of the circumstellar environment. In contrast, the 
optical properties in regions of low absorption may be crucial 
for grain temperatures and radiative acceleration, if they coin- 
cide with regions of high stellar flux. 

As an example, consider the optical data for AI2O3 pro- 
vided by Koike et aT] ( | 1995[ > and the extrapolated data used by 
Woitke (2006). The slopes of the absorption coefficients in the 



near-infrared differ dramatically for these two sets of data, with 
p ~ 1.7 and p as -2.1, respectively. The resulting grain tem- 
peratures at a given distance from the star will therefore be very 
different, affecting the condensation distance. The numerical re- 
sults in Woitke ( 2006 ) suggests that AI2O3 grains form as close 
as 1 .5 R„ to the stellar surface, whereas a simple estimate for the 



data from Koike et al. ( 1995}, using Eq. ( 10 1, results in a conden- 
sation distance of about 3 R* for T„ = 2500 K. However, setting 
aside uncertainties in optical data and the condensation distance, 
AI2O3 is an unlikely candidate since the elemental abundance of 
Al is 1-2 orders of magnitude lower than for most of the limiting 
elements in Tab. Q] 

It is likely that several different kinds of dust species will 
form close to the stellar photosphere and not all of them will 
fulfill the criteria concerning abundance and cross-section that 
are required to trigger outflows. However, some of them may 
act as seed particles for the actual wind-drivers. We might also 
expect silicate grains to be enriched by iron further out in the 
wind, even if such grains cannot form close to the stellar surface. 

To summarize the above discussion we compile the criteria 
for wind-drivers in Tab. [3] The criterion concerning the cross- 
section has been divided into two columns, one for the small 
particle limit (SPL) and one for bigger grains where scattering 
may play a role (BG). The column with cross-sections for grain 
sizes beyond the small particle limit has been filled in when the 
contribution by scattering may be the decisive factor for mak- 
ing a material a potential wind-driver (as in the case of SiCh, 
MgSiC>3 and M^SiCu). The chemical and optical properties of 
AI2O3, TiC<2 and MgSi03 place them just outside the part of the 
parameter space producing outflows (see Fig.[TJand|4|. However, 
they have been listed as satisfying the condition of condensation 
distance, since the boundary between the wind/no-wind zones 
might be slightly shifted due to varying pulsation amplitudes 
and stellar parameters (see Fig. [TJ. More examples of p and T c 
values for dust species that are not considered potential wind- 
drivers, due to low abundances or large condensation distances, 
are given in Tab. [4] 

We have included amorphous carbon in Tab.[T]and[3]for com- 
parison because it is a known wind-driver for C-type AGB stars. 
Amorphous carbon (amC) is very opaque in the relevant wave- 
length region, it has a high condensation temperature and the 



Table 4. More examples of p and T c values (column 2 and 3) 
for dust species in the circumstellar environments of AGB stars, 
as inferred from observations (Mo lster et al.[ 20 1 0) > . In column 4 
the condensation distances R c , using Eq. ( 10 1, are listed. 



Material 


P 


T c [K] 


RJR, 


Ref. (r c and p) 


MgAl 2 4 


-1.2 


1150 


1.7 


1, 3. 


Mgo. 5 Feo.5Si0 3 


2.4 


1100 


9.9 


1, 2 


Fe 2 Si0 4 


1.7 


1050 


8.2 


1. 3 


Mgo.5Feo.5O 


1.9 


1000 


10.4 


1,2. 


FeO 


1.6 


900 


12.0 


1,2. 


SiC 


0.3 


1100 


3.7 


1, 4. C-rich atmos. 



References. (1) Gail 2010); (2) |Henning etai1fT9"9"5) ; (3) |Zeidler et aT] 



p0TT] >; (4)|Pegourie 



1988 



limiting element is abundant in carbon-rich atmospheres. There 
is, however, practically no carbon available for forming dust par- 
ticles in M-type AGB stars, since most carbon atoms are locked 
in the tightly bound CO-molecule. The optical and chemical 
properties of this dust species (p « 1,2 and T c * 1700 K) place it 
in a different region of the parameter space than the grain materi- 
als discussed above. Therefore amC grains will affect the spectra 
in a different way than e.g. Fe-free silicate particles, as will be 
discussed in a forthcoming paper (Bladh et al. 2012[ ). 



7. Summary and conclusions 

Observations of molecular lines in the inner part of the atmo- 
sphere provide constraints on how far shock waves can levitate 
gas (i.e. a few stellar radii). A necessary condition for the occur- 
rence of dust-driven outflows is that the condensation distance 
falls within this range. However, it is not sufficient that grains of 
a certain material can form within reach of the levitated gas in 
order to trigger a wind (criterion (a) in Sec 2.3 ). The individual 



dust species also need radiative cross-sections and abundances 
such that they can provide the radiative acceleration necessary 
to cause a general outflow (criteria (b) and (c) in Sec 2.3 1. 

To investigate how a close stellar environment affects dust 
formation we have a constructed a grid of detailed atmosphere 
and wind models, including a frequency-dependent treatment of 
radiative transfer for the gas and dust components. We see dis- 
tinct trends in the numerical results concerning what combina- 
tions of optical and chemical dust properties (here quantified by 
the parameters T c and p) are required to trigger an outflow. Dust 
species with a low condensation temperature and a near-infrared 
absorption coefficient that decreases strongly with wavelength 
will not condense close enough to the stellar surface to be con- 
sidered as potential wind-drivers. The grain temperature, and 
consequently the condensation distance, is strongly influenced 
by the wavelength dependence of the absorption coefficient. A 
simple estimate (see Eq. ( [T0| )) of how the interaction between a 
non-grey dust opacity and the stellar radiation field affects the 
condensation distance shows good qualitative agreement with 
the detailed numerical models. However, the quantitive results 
differ due to the inherent physical simplifications in the analyti- 
cal estimate. 

A frequency-dependent treatment of the dust opacity is abso- 
lutely necessary when calculating dust yields from AGB stars. A 
grey treatment of the dust opacity (corresponding to p = 0) will 
underestimate the condensation distance for dust species with a 
power law coefficient p > while the opposite is true for dust 
species where p < 0, producing misleading results regarding the 
overall dust composition. 
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Table 3. Compilation of properties for a few selected dust species. 



Grain material 


Condensation 


Abundance (b) 


Cross-section (c) 


Wind-driver 




distance (a) 




SPL 


BG 




Fe 


NO 


YES 


YES 




NO 


A1 2 3 


YES 


NO 


NO 




NO 


Ti0 2 


YES 


NO 


NO 




NO 


Si0 2 


? 


YES 


NO 


? 


? 


MgSiOa 


YES 


YES 


NO 


7 


? 


Mg,Si0 4 


YES 


YES 


NO 


YES 


YES 


MgFeSi0 4 


NO 


YES 


YES 




NO 


amC 


YES 


YES 


YES 


YES 


YES 



Notes. A compilation of the criteria mentioned in Sec [53] for the dust species listed in Tab. [T] The question mark in condensation distance and 
cross-section for Si02 and AI2O3 is due to uncertainties in the optical data. Note that amC is a wind-driver in C-type AGB stars whereas all other 
dust species are considered in the context of M-type AGB stars. For figures showing dust opacities beyond the small particle limit, see [Hofner 
(2008! and |Mattsson & H6fner| < |2011} for Mg 2 Si0 4 and amC, respectively. 



As for specific grain materials, the wavelength dependence 
of the near-infrared absorption coefficients for metallic iron and 
Fe-bearing silicates, together with their low condensation tem- 
peratures, prevent them from forming close to the stellar surface, 
in accordance with earlier results by Woitke (2006 ). This means 
that neither metallic iron or Fe-bearing silicates are viable wind- 
drivers. T1O2 is also excluded as a candidate, but in this case it is 
a consequence of the low abundance of the limiting element, Ti, 
and the low cross-section. Other cases, such a SiC»2 and AI2O3, 
are less clear-cut due to uncertainties in the dust data and further 
work is needed. In particular, precise optical data in the near- 
infrared where AGB stars emit most of their flux is required for 
a reliable determination of grain temperatures. 

A dust species that has been shown to work as a wind-driver 
if the grains grow to sizes around 0.1-1 micron is Mg2SiC>4 
( Hofnerj |2008 ). It consists of relatively abundant elements and 
the condensation temperature together with the wavelength- 
dependent dust opacity allows it to form close to the stellar 
surface (~ 2R t ). These grains are very transparent in the near- 
infrared in the small particle limit, but if the grains reach sizes 
comparable to the wavelength of the stellar flux maximum and 
scattering is taken into account, the overall cross-section can 
provide enough radiative acceleration to trigger a wind. This sce- 
nario is supported by recent observations of grains in the critical 
size-range (i.e. ~ 0.3 /urn) in the close stellar environment ( |Norris| 
[eTaTj[20T2l i. 
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Appendix A: Radiative acceleration of the wind 

The collective opacity for dust particles of radius a gl is given by 



P ° P 



7T QacJA, flgr) a 

a„ r n 



gr"gr> 



(A.l) 



where « gr is the number density of grains in the atmosphere, Q dcc 
is the efficiency and p is the total mass density. The factor aLrip 
is related to the fraction of atmospheric volume occupied by dust 
particles. This quantity can be formulated in terms of the volume 
of a monomer, the basic building block of the grain material, 
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Here we have expressed the volume of the monomer V mon in 
terms of the atomic weight of the monomer A mon , the proton 
mass m p and the bulk density p gl . The number of monomers con- 
tained in all grains in a volume of atmosphere n mon can be ex- 
pressed by the abundance of the limiting element eij m , the num- 
ber of corresponding atoms in the monomer s, the fraction of the 
limiting element condensed into dust particles f c and the total 
number density of H atoms « H . Using « H = p/(l + 4e He ), we 
obtain Eq. |6]) for the dust opacity 
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Setting Pd = Pc and solving for R c — r(T c ) results in the expres- 
sion for the condensation distance seen in Eq. dlOj, 
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For more details, see |Lamers & Cassinelli| (l999) 

Appendix C: Levitation distance 

Conservation of energy, with kinetic energy fully converted into 
potential energy, for a fluid element of mass m, with initial ve- 
locity m ( ) at a distance Rq from the star, results in the following 
expression for the levitation distance R mdx , 
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where the escape velocity at the stellar surface is given by 
Mesc = (2M t G/R t ) l/2 . 



Appendix B: Condensation distance 

The grain temperature T$ is determined by the condition of ra- 



diative equilibrium (see Sect. 3.1 

fabs ; J-/-fabs,s5(Pd) = 0. 



(B.l) 



Assuming an optically thin atmosphere, where the incident in- 
tensity on the grains is direct star light, the mean intensity J A can 
be approximated by a geometrically diluted Planck distribution, 

J A = W(r)B A (T t ). (B.2) 

In this expression W(r) is the geometric dilution factor 

W(r)= ^l-(RJr) 2 ], (B.3) 

which reduces to W(r) = (R t /2r) 2 for r » R t . The condition of 
radiative equilibrium can thus be reformulated as 



(B.4) 



where the source function S (Pa) is given by B A (Ti). Assuming 
r » R* and solving for the grain temperature (Pd) gives 
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If the absorption coefficient can be approximated with a power 
law function, /c a b s ~ X~ p , the factor with the Planck mean opaci- 
ties can be simplified accordingly, 



'Kabs,B(P»)\ Pi' 



\Kabs,B(Td)l Pj' 

resulting in the following grain temperature distribution 
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